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Abstract 



We show that the assumption of type II seesaw mechanism for small neutrino 
masses coupled with b — r mass unification in a minimal SUSY SO(10) model 
leads not only to a natural understanding of large atmospheric mixing angle 
(#23) among neutrinos, as recently pointed out, but also to large solar angle 
(#12) and a small #13 = U e 3 as required to fit observations. This is therefore a 
minimal, completely realistic grand unified model for all low energy observa- 
tions that naturally explains the diverse mixing patterns between the quark 
and leptons without any additional inputs such as extra global symmetries. 
The proposed long baseline neutrino experiments will provide a crucial test 
of this model since it predicts U e s ~ 0.16 for the allowed range of parameters. 



1 



I. INTRODUCTION 



The various neutrino oscillation experiments such as those involving solar and atmo- 
spheric neutrinos as well as the KEK and KAMLAND experiments that involve laboratory 
produced neutrinos have now produced quite convincing evidence that neutrinos have mass 
and they mix among themselves. Although the neutrinos have a great deal of similarity with 
quarks as far as the weak interactions go, the oscillation results have revealed a profound 
difference i.e. two of the three neutrino mixings are very large whereas all quark mixing 
angles are very small. Understanding this difference is a major challenge for theoretical 
particle physics today. The problem becomes particularly acute in models that unify quarks 
and leptons such as the SO(10) grand unified models [1], which are considered as prime 
candidates for describing neutrino masses. 

Some of the reasons that make SO(10) models so attractive as grand unification theories 
of nature are the following: (i) in SO(10) model, all fermions can be part of a single spinor 
representation; (ii) it contains the left-right symmetric unification group SU(2) L x SU (2) R x 
U(1)b-l x SU(3) c [2] which provides a more satisfactory way to understand the origin of 
parity violation in Nature and (iii) finally, perhaps the single most important reason is 
the natural understanding of small neutrino masses via the seesaw mechanism [3] since the 
single spinor representation discussed above that contains all the standard model fermions 
also contains the right handed neutrino, needed in implementing the seesaw mechanism. 

A closer look at the details of neutrino oscillation data in fact provides one more com- 
pelling reason for the SO(10) model: in order to understand the atmospheric neutrino data, 
we need the heaviest neutrino mass to be larger than 0.05 eV. The seesaw formula i.e. 
m v ~ ^ (where mo is the neutrino Dirac mass and M R is the mass of the right handed 
neutrino contributing to the neutrino mass needed to understand the atmospheric data) then 
tells us that there must be one right handed M R < 10 15 GeV. This value is considerably 
smaller than the Planck mass and therefore one is faced with a new hierarchy problem sim- 
ilar to the corresponding problem of the standard model. However, it was pointed out long 
ago [4] that the Majorana mass of the RH neutrino owes its origin to the breaking of local 
B-L symmetry which implies that M R ~ M B _ L . Local B-L symmetry therefore provides a 
natural way to understand the smallness of the RH neutrino mass compared to M P£ . What 
is very interesting is that SO(10) group also contains the local B-L as a subgroup. 

SO (10) model, despite its attractiveness for understanding the overall scale of neutrino 
masses, runs into a potential trouble in providing an understanding of the observed mixings. 
The problem arises from the fact that SO(10) also contains the quark-lepton unification 
577(4) c group of Pati and Salam, which in the simplest approximation leads to equal quark 
and lepton mixing angles and one needs to make further assumptions to get a handle on the 
mixings [5]. An obvious conceptual problem is that if one of these models is ruled out by 
data, one would not be able to tell whether it is the SO(10) unification which is "at fault" 
or it is one of the assumptions used to derive neutrino mixings. 

A different approach to this issue was taken in ref. [6]. The idea was to avoid the use 
of any symmetries beyond the gauge symmetry, in this case SO (10) and use the minimal 
set of Higgs fields that can break the group down to the standard model and give mass 
to the fermions. It turns out that if we choose Higgs fields in 10 and 126, the Yukawa 
superpotential contains enough parameters to fit the observed fermion masses and mixings 
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of the standard model. It was observed in ref. [6], that in this model the neutrino masses and 
mixings are completely predicted upto an overall scale, when one uses the seesaw mechanism 
which is part of the SO(10) model 1 . To break the gauge symmetry fully, two additional Higgs 
multiplets belonging to 45+54 are included. They do not contribute to the fermion masses, 
leaving the conclusions on neutrino masses unchanged. 

An additional appeal of breaking B-L symmetry of SUSY SO(10) by an 126, as opposed 
to by 16 Higgs, is that it automatically leaves R-parity as an exact aymmetry and thereby 
explains, why neutralino is a stable dark matter [7,8]. This is because the submultiplet of 
126 that breaks B-L carries B-L =2. Therefore R-parity (defined by Rp = (-lf( B - L )+ 2S ) 
quantum number of this field is even and therefore, its vev leaves R-parity unbroken. In 
contrast in models where B-L is broken by a 16-plet of Higgs, the B-L symmetry is broken by 
one unit and without any additional symmetries (e.g. matter parity), neutralino is unstable 
and cannot therefore serve as a dark matter. Of course, if a fundamental theory e.g. a 
superstring theory that led to an SO(10) model with appropriate additional symmetries 
that guarantee the stability of neutralino was known, then the above objection to an 16 
Higgs would not apply. 

The initial analyses of neutrino mixings within the framework of Ref. [6] used the simple 
seesaw formula (type I seesaw) and are now in disagreement with data. In subsequent papers 
[9-13], this idea has been analysed ( in some cases by including more than one 10 Higgses) 
to see how close one can come close to the observed neutrino parameters. The conclusion 
now appears to be that one needs CP violating phases to achieve this goal, as noted in [13]. 

A way out of this problem is to use the type II seesaw mechanism, as was initially done 
in [11], where an induced triplet vev is added to the usual type I seesaw formula arising from 
the RH neutrino intermediate state. In models which have asymptotic parity symmetry such 
as left-right or SO(10) models, type II seesaw arises if both parity and B-L symmetry are 
broken at the same scale. 

A very interesting point about this approach has been noted in a recent paper [14], 
where it has been shown that if we restrict ourselves to the 2-3 sector of the model and use 
the type II seesaw mechanism, then the b — r unification of supersymmetric grand unified 
theories leads to a neutrino Majorana mass matrix which explains the large — v T mixing 
angle needed to understand atmospheric neutrino data. The important point is that no 
symmetries are needed to get this result. 

To see how this result is very generic to SO(10) models with 126 vev, note that: (i) in 
SO (10) model, the neutrino mass matrix is given by the type II seesaw formula, with two 
contributions: one coming from the righthanded neutrino intermediate state and another 
coming from an induced triplet vev generic to these models, as already mentioned; and (ii) 
that for certain range of parameters, the induced triplet vev term can dominate the neutrino 
mass matrix. As was shown in [11], under these circumstances, one gets a sumrule 

M v = a(M e -M d ) (1) 



1 This is to be contrasted with the SU(5) case where the minimal Higgs set needed to break 
the gauge symmetry i.e. 5+24 Higgses lead to the mass relation m e /m^ = md/m s that is in 
contradiction with observations. 
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In [14], it was observed that since this relation is valid at the seesaw scale, one must use the 
extrapolated quark and lepton masses in the formula. The fact that at or near the GUT 
scale m b /m T ~ 1 — 1.2 depending on the value of tan/3, implies that the 3-3 element of 
the M. v which is proportional to m b — m T is of the same order as the off diagonal elements 
of the Mnu in the 2-3 subsector leading top the largeness of the atmospheric mixing angle 
without any further assumptions. 

It is however essential to do a complete three generation analysis of this model if this 
important observation is to lead to a realistic SO(10) model for understanding all neutrino 
mixings. In fact, since the model has no free parameters, it is a priori not obvious that 
within this framework one would simultaneously get a large solar mixing angle and a small 
9 13 = U e s as well as the correct value for the ratio Arn^/ Am\. It is the goal of this paper 
to analyze this question. 

We find that for a narrow range of the quark masses, the model does indeed lead to 
the correct mass difference-squares as well as large mixing angles 9 12 and 6*23 needed to 
understand the solar and atmospheric neutrino oscillations while at the same time keeping 
U e 3 < 0.16 as required by the reactor data. We find that for this entire range, the reactor 
angle C/ e 3 — 0.16 providing a clear way to test the model in the proposed long base line 
neutrino experiments. 

This paper is organized as follows: in sec. II, we review the basic equations of the model 
for the supersymmetric case; in sec. Ill, we explain the conditions under which the induced 
triplet vev dominates neutrino masses; in sec. IV, we describe our method for solving the 
equations to predict the neutrino parameters. In sec. V, we give our results for the neutrino 
masses and the mixing angles. In sec. VI, we give our conclusions. In this paper, we ignore 
the CP violating phases. 

II. THE MASS SUMRULES FOR MINIMAL SO(10) 

We consider supersymmetric SO(10) group with the Higgs fields belonging to the rep- 
resentations 45+54 for breaking SO(10) group down to the left-right symmetric group 
SU(2) L x SU(2) R x U{l) B -L x SU(3) c and the minimal Higgs set 10+126+126 that cou- 
ple to matter and also break the SU(2) L x SU(2) R x U(1) B -l x SU(3) c group down to 
SU(3) C x U(l) em . It is the latter set i.e. 10 © 126 which is crucial to our discussion 
of fermion masses. The first stage of the symmetry breaking therefore could have been 
accomplished by alternative Higgs multiplets e.g. by 210 of Higgs without effecting our 
results. As has been noted earlier [6,10], the set 10+126 which couple to matter contain 
two pairs of MSSM Higgs doublets belonging to (2,2,1) and (2,2,15) submultiplets (under 
SU(2) L x SU(2) R x SU{A) C subgroup of SO(10)). We denote the two pairs by <f) u4 and A u4 . 
At the GUT scale, by some doublet-triplet splitting mechanism these two pairs reduce to 
the MSSM Higgs pair (H u , H d ), which can be expressed in terms of the <fi and A as follows: 

H u = cosa u <p u + sina u A u (2) 
H d = cosa d (f) d + sina d A d 

The details of the doublet-triplet splitting mechanism that leads to the above equation are 
not relevant for what follows and we do not discuss it here. As in the case of MSSM, we will 
assume that the Higgs doublets H U:d have the vevs < H® >= vsinfi and < >= vcosf3. 



4 



In orders to discuss fermion masses in this model, we start with the SO(10) invariant 
superpotential giving the Yukawa couplings of the 16 dimensional matter spinor ipi (where 
i,j denote generations) with the Higgs fields H w = 10 and A = 126. 

W Y = hijipitpjHn) + fij^itpjA (3) 

SO(10) invariance implies that h and / are symmetric matrices. We ignore the small effects 
coming from the higher dimensional operators. Below the B-L breaking (seesaw) scale, we 
can write the superpotential terms for the charged fermion Yukawa couplings as: 

W = h u QH u u c + h d QH d d c + h e LH d e c + fiH u H d (4) 

where 

h u = hcosa u + fsina u (5) 
h d = hcosa d + fsina d 
h e = hcosa d — 3fsina d 

In general a u ^ a d and this difference is responsible for nonzero CKM mixing angles. In 
terms of the GUT scale Yukawa couplings, one can write the fermion mass matrices at the 
seesaw scale as: 

M u = (h + f) (6) 
M d = {hri + fr 2 ) 
M e = (hn-3r 2 f) 
M v u = (h-3f) 

where 

h = hcosa u sinf3 (7) 

/ = f sina u sinf3 
cosa d 

T\ = COtp 

COSCi u 

sina d 

r 2 = — COtp 

sina u 

To count the number of parameters describing the fermion sector, we first ignore CP phases 
and choose a basis where h is diagonal. Since / is symmetric, we have a total of nine 
parameters from the couplings and including a u>d and (5 gives us a total of twelve parameters. 
All these parameters can be determined by fitting the the six quark masses, three lepton 
masses and three CKM angles. This enables a complete determination of the neutrino masses 
upto an overall scale related to the B-L symmetry breaking and the three mixing angles. 
The model is therefore completely predictive in the neutrino sector. 

In order to determine the neutrino masses and mixings, one uses the seesaw mechanism. 
As noted in the introduction, most previous works on this model except the works in Ref. 
[11,14] used the type one seesaw mechanism where the neutrino mass matrix is given by the 
formula: 



5 



M v = -M v uM^M T v u 



(8) 



where M^ R = fvs-L- On the other hand it is well known that in asymptotically parity 
conserving theories including SO(10), the true seesaw formula [15] (called type II here and 
in literature) has a second term which arises from an induced SU{2)l triplet vev or from 
higher dimensional terms involving left doublets: 

M v = ]a L -M vD M N \M T vD (9) 

where ol = A v2 , where v is the SU{2) L breaking scale and A is a combination of 
parameters in the Higgs potential. The type II seesaw formula has two new parameters 
and vb-l instead of one in the type I case; but for different ranges of initial parameters 
in Higgs potential, either the first or the second term can be made to dominate [16]. We 
will work within the assumption that it is the first term that dominates the seesaw formula 
[14]. This situation can arise when ^> v 2 /vb-l (or A which is a ratio of scalar coupling 
parameters is much larger than one). We elaborate the circumstances when this happens in 
sec. III. 

It was noted in several papers [6,10-12] that if one uses type I seesaw, this model cannot 
produce two large neutrino mixing angles if CP phases are ignored. It has been shown 
recently [13] that once the CP phases are included, one can get bi-large neutrino mixing 
pattern. But here we search for two large mixing angle solutions without invoking the 
CP phase. As noted in [14], if one uses type II seesaw and assumes further that the first 
term dominates, large atmospheric mixing angle follows naturally as a consequence of b-r 
unification. Whether this also simultaneously yields a large solar mixing angle and a small 
U e 3 remained an open question. In this present paper we present a detailed analysis of this 
idea in a full three generation model and show that for a very narrow range of quark masses, 
this model also predicts large solar mixing angle as well as the correct solar mass splitting 
together with a small U e ^. Since the model has hierarchical mass pattern, the atmospheric 
mass difference-squared is related to the overall scale of the triplet vev in the type II seesaw 
formula and cannot be predicted. 

III. DOMINANCE OF INDUCED TRIPLET VEV 

In this section, we would like to discuss the parameter range where the induced triplet 
vev term dominates the neutrino mass matrix. To elucidate this let us discuss the origin of 
the triplet vev in our minimal SO(10) model. First we note the decomposition of the 126 
under the group SU(2) L x SU{2) R x SU{A) C : 

T26 = (1,1, 6) ©(2, 2, 15) ©(3, 1,10) 0(1, 3, 10) (10) 

The SU{2)l triplet that contributes to the type II seesaw formula is contained in the mul- 
tiplet Al = (3,1,10) and it couples to the left handed multiplet ip = (2,1,4) of the 16 
dimensional SO(10) spinor that contains the matter fermions i.e. ^l^l^-l- On the other 
hand the mass of the RH neutrinos comes from the coupling of Ar = (1,3, 10) submultiplet 
of 126 to the right handed fermion multiplet ip c L = (1,2,4) i.e. i\) c L i\) c L A R . 
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The vev of the neutral member of A# breaks the B-L symmetry and gives mass to the RH 
neutrinos and generates the second term in the type II seesaw formula. To see how the A° L vev 
arises, note that the general superpotential of the model contains terms of type Ail26 2 • 54 
and A 2 10 • 10 • 54. In the Higgs potential, this generates a term (from I-F54I 2 ) of the form 
10 • 10 • 126 • 126. In this term, there is a term of the form 0(2, 2, 1) 2 A L (3, 1, 10)A fl (l, 3, 10) 
with a coefficient A1A2. Furthermore, in the Higgs potential, there is a mass term for 
Ax, (3, 1, 10) of the form n^ + X^vlf, where vu is the GUT scale. On minimizing the potential, 

these two terms lead to a vev for the SU(2) L triplet a L =< A° >~ Xl ^ V +\^ L ■ ^ 1S now 
clear that if we choose A3 such that fi^ + XsVu <C v%_ L , then the entries in the second matrix 
in the type II seesaw formula can much smaller than and Eq.l holds. We will work in 
the domain of the parameter space where this happens. 

If the triplet vev contribution to the neutrino mass matrix dominates in the type II 
seesaw formula, then Equation (6) can be used to derive the sumrule 

M u = a(M e -M d ) (11) 

Using this equation in second and third generation sector, one can understand the results of 
[14] in a heuristic manner as follows. The known hierarchical structure of quark and lepton 
masses as well as the known small mixings for quarks suggest that have the following 
pattern 

M e * m T ( 6 ^ j) (12) 

and 

M d * m b ( td > 1 6 f) (13) 

where C 1. It is then clear that if there is approximate b — r unification as it appears 
to be the case if the theory below the B-L breaking scale is MSSM, then in — M d 
matrix, there is a high degree of cancellation in the 33 entry making this entry comparable 
to all the other entries in this matrix. The atmospheric mixing angle which is given by 
tan20A — (m T e< i2 — rn h e d ^)l{rn h — m T + m b €d t i — mie^i) becomes very large at the B-L 
breaking scale. Since renormalization played an important role in obtaining this result, one 
must ask what happens to the neutrino mixings once they are extrapolated to the weak 
scale [17]. It is well known [17] that for the case of normal hierarchy for neutrino masses 
as is the case here, the MSSM RGE's do not change the mixing angles very much and the 
seesaw scale result persists at the weak scale with only minor changes. 

IV. DETAILS OF CALCULATION 

In this section, we outline our method for determining the neutrino mixing parameters. 
For this purpose, we first note that the matrices h and / in Eq. (6) can be eliminated in 
terms the mass matrices M u>d so that we have a sumrule involving the three mass matrices 
M Ut d,e- Before giving the sum rule, we note that we will work in a basis where M d is diagonal 
and is given by M u = V T ■ • V (where is the diagonal mass matrix of up type quark 
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and V is the CKM mixing matrix). This can be done without any loss of generality. We 
also introduce a new set of matrices Mi jU> d where M = m 3 being the third family mass 
for the corresponding flavor. The sumrule for charged lepton matrices is given by: 

kMi = rM d + M u 

where k and r are functions of 7*1,2 (given in sec. II) and fermion masses as follows: 

r 2 - 7*1 m T 



k 



r 2 + 3ri rn b 
4r 1 r 2 m t 



(14) 

(15) 
(16) 



M v = ai—Ma-W) 
m-r 



(17) 



These relations are valid at the B-L breaking scale vb-l- The advantage of working with 
M rather tha M is that the 33 elements of all M^ u ^ matrices are either one or of order 
one; so we expect solutions for k and r also of order one. Furthermore since the formula 
for M. v involves only Mi and Md, b — r unification helps to see the cancellation in the 33 
element of hA v somewhat more easily. At the same time the 23 element of A4 U receives 
only one contribution from since in our basis Md is diagonal. These two results lead to 
atmospheric mixing angle being large [14]. 

To carry out the calculations, we have to solve for the two unknowns k and r using the 
low energy inputs from the quark and charged lepton sectors. To obtain a perturbative 
estimate of these parameters, we decompose rM d + M u as: 



(18) 



where e i: a <C 1 and y. In this analytical approach, our procedure will be to find the 

eigenvalues of (18) by perturbation method and match them to the known leptonic masses at 
the B-L scale. The advantage of this decomposition is that it allows a nice perturbative de- 
termination of the eigenvalues analytically without having to resort to immediate numerical 
analysis. We will compare our results with the numerical evaluation using Mathematica. 
The i th eigenvalue \ = A- - 1 + A- 2 - 1 is found to be 

= x (19) 
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V? + 4e| 



A 
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+ z + yj(z-y)* + 4e 



2 + ^ + ^0(10" 



{Z61 - ae 2j 



0(1(T 



X (2) 



Ai 



(2) 



~O(10~ 2 )Af 
-(Af 



(2) + A? } ) 
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We consider only cases where y ~ 10~ 2 and z > 0.1. Within this regime, the unperturbed 
2nd and 3rd lepton masses are accurate up to a few %. However, the higher order electron 
mass correction is big and so the perturbation formula breaks down for this case. We 
therefore use the perturbation technique for the second and third generation masses but 
use the determinant to find that for the first generation. As mentioned, we will check the 
validity of perturbation result using numerical methods. 

Taking determinant of the above equation 18, we find that the three charged lepton 
masses are related as follows: 

/c 3 m e m M = xyz — rre 2 , — ya 2 — zt\ + 2ae^ 2 (20) 

km^ = A 2 ~ A 2 0) (21) 



k = A 3 ~ A 3 0) ~ z + ^ (22) 

z 



We now solve the above equation by substituting x,y,z, a,ei,e 2 with the corresponding ele- 
ments in the matrix rMd + M u . From eq.(20), and find 

fc(l + m M ) = y + z (23) 

k = z + — (24) 

z 

Since Eq.18 tells us that z = 1 + r and y = rs + c, we can use the above two equations 
to determine the parameters k and r, which we can then use to find neutrino masses and 
mixings. We find k and r to be 



(s + c - 2m^) ± y/(s- c) 2 - 4(m^ - s) (1 + m^)t\ 
2(m M - s) 

(l + s)r + l + c 



(25) 



k = 



1 + m M 



and a consistency relation for the d-quark mass 



= k^fh e fh^ + ze 2 + ya 2 - 2ae 1 e 2 - u(yz - e 2 .) 
r(^ - e|) 

In order to get a rough feeling for the way the maximal neutrino mixings arise, let us 
diagonahze the charged lepton mass matrix given in Eq. 14 and write the neutrino mass 
matrix in this basis: 



M v = a(^-U}M d U e -Mt D ) (27) 

m. 



Where Mf is the diagonal charged lepton mass matrix with r mass is 1. C/j is the rotation 
matrix diagonahze charged lepton mass. U\ can be written approximately as 



U, 



( 1 5 l 5 2 \ 

Ai cos0 sin0 
\ A 2 — sin <j) cos J 



(28) 
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where 

tan^A^,. (29) 
z — \ 2 

The parameters Si and Aj are given to lowest order in perturbation theory by 

s _ 61 cos - a sin 
kfn^ — x 
€\ sin + a cos 
k — x 

Ai = —Si cos — S2 sin 
A 2 = 5i sin — 5 2 cos 

Using these parameters and neglecting small terms due to Si and S 2 multiplying light quark 
masses, we find that 



M v ~ 



( m d — m e + m s A 2 + m h A\ m s Ai cos — m b A 2 sin m s A x sin + m b A 2 cos ^ 

m s Ai cos0 — m;,A 2 sin m s — + mj sin 2 — m;,sin0 

y m s Ai sin0 + m h A 2 cos0 — m b sin0 — m fe sin 2 + m b — m T y 



(31) 

We now find the following analytic expression for the atmospheric mixing angle from Eq. 
31 to leading order ignoring small terms to be: 

^^^^r—T^T^ ( 32 ) 



Q = 



q + V<? 2 + 4 
2mt sin 2 + (m T — mi,) + (m s — m M ) 
to& sin 



For |g| < 1, we get sin 2 29a > 0.8. We see that 6 — r unification i.e. m b ~ m r and 
mf,sin0 ~ (rrib — m T ) are important to get a large 9a- Also we need to have m s < and 



> 0. 



V. PREDICTIONS FOR NEUTRINO MASSES AND MIXINGS 

In order to obtain the predictions for neutrino masses and mixings in our model, we will 
need the values of quark masses and mixings at the seesaw scale. Experiments determine 
these input parameters near the GeV scale and they need to be extrapolated to the B-L 
scale which is near 10 15 -10 16 GeV where our Equation (6) is valid. Taking the values for the 
quark masses and mixings at the B-L scale we can determine k and r approximately. We 
will use this determination of k and r to solve for neutrino masses and mixings using the 
relation in Eq.17. We will also compare our results with a direct numerical scan of the Eq. 
14 i.e. not using perturbation method to obtain k and r. Results obtained by both methods 
are in agreement. 
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In our model, the theory below the B-L breaking scale is the MSSM whose effect on 
fermion mass extrapolation is a well studied problem [18]. We will use the two loop analysis 
in the paper by Das and Parida [18] in our analysis. Our strategy will be to take the values 
of the quark masses at the scale v B -l — 10 16 GeV given in [18]. In Table I, we give the 
input values of masses and mixings for values of the MSSM parameter tan(3 = 10 and 55. 



input observable 


tanp = 10 


tanfi = 55 


m u (MeV) 


79 +u - i3 


72 +u - 12 

u - ' z -0.14 


m c (MeV) 


210.32™ 


210.50+^-^ 


m t (GeV) 


82.43lSre 


95.14™ 


m d (MeV) 


l.OU_ 23 


1-49^ 


m s (MeV) 


29.94+^ 


29.81^ 


m b (GeV) 


1 nfi +u - 14 


1 41+U.48 


m e (MeV) 


0.3585 


0.3565 


m M (MeV) 


75 671 ci+u- 05 ^ 

' °- u ' lo -0.0501 


75.2938i^ 15 


m T (GeV) 


i o Q 99+U.UU13 

i.zyzz_ .ooi2 


1 fi oqo+U.U443 
l.OZyz_ 0294 



Table I: The extrapolated values of quark and lepton masses at the GUT scale from the 
last reference in [18]. We have kept the errors to only two significant figures in the quark 
masses. 

For the mixing angles at GUT scale, we take: 

/ 0.974836 0.222899 -0.00319129 \ 
Vckm = -0.222638 0.974217 0.0365224 (33) 
V 0.0112498 -0.0348928 0.999328 / 

In the first perturbative method, we use the above input values to obtain k and r using Eq. 
25 and search for values around them that give a good fit to charged lepton masses and 
then use them in Eq.17 to derive the neutrino masses and the three mixing angles: sin 2 29 Q , 
sin 2 26 A and U e3 . The best fit range for k,r are -.78 < r < -.74 and 0.23 < k < .26. 
We also do a direct numerical solution. Both the results are in agreement. (We ignore CP 
violation in this work.) 

Note that the sign of a fermion is not physical, which leads to several choices for the sign 
of fermion masses that we have put into our search for solutions. The solutions we present 
here correspond to m e ^ T>b j > and m c ^ s < upto an overall sign. 

Our results are displayed in Fig. 1-3 for the case of the supersymmetry parameter 
tanfi = 10. In these figures, we have restricted ourselves to the range of quark masses for 
which the atmospheric mixing angle svtPIQa ^ 0.8. (For presently preferred range of values 
of sin 2 29A from experiments, see [19]). We then present the predictions for sin 2 29 Q , Am|, 
and U e3 for the allowed range sin 2 29 A in Fig.l, 2 and 3 respectively. The spread in the 
predictions come from uncertainties in the s, c and the 6-quark masses. Note two important 
predictions: (i) sin 2 26 & > 0.91 and U e3 ~ ±0.16. The present allowed range for the solar 
mixing angle is 0.7 < sin 2 29 Q < 0.99 at 3a level [19,20]. The solutions for the neutrino 
mixing angles are sensitive to the b quark mass. 
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It is important to note that this model predicts the U e s value very close to the present 
experimentally allowed upper limit and can therefore be tested in the planned long base line 
experiments which are expected to probe U e3 down to the level of ~ 0.05 [21,22]. Our model 
would also prefer a value of the sin 2 29 A below 0.9, which can also be used to test the model. 
For instance, the JHF-Kamioka neutrino experiment [22] is projecting a possible accuracy in 
the measurement of sin 2 20A down to the level of 0.01 and can provide a test of this model. 




0.8 0.82 0.84 0.86 0.88 0.9 



Sin 2 2e A 

FIG. 1. The figure shows the predictions for sin 2 26 & and sin 2 20A for the range of quark masses 
in table I. Note that sin 2 20 Q > 0.9 and sin 2 29 A < 0.9 



As tan/3 increases, the allowed values for the neutrino mixings and masses fall into an 
even narrower range. 

VI. CONCLUSION 

In summary, we find that a minimal SO(10) model with single 10 and 126 Higgs coupling 
to matter is a completely predictive model for neutrino masses and can provide an excellent 
description of the presently favored patterns for neutrino masses and mixings required by 
data. The only assumption needed is that the SU{2)l triplet vev dominates the neutrino 
masses. No global symmetries are invoked to generate the neutrino mass pattern, unlike 
most models that employ the 16 Higgs to break the B-L symmetry. The model predicts 
a hierarchical mass pattern for neutrinos and a value of U e 3 ~ 0.16, both of which can be 
tested in upcoming long baseline neutrino experiments. The atmospheric mixing angle is 
found to be between 0.8 and 0.9 which is also a testable prediction of the model. In our 
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0.8 0.82 0.84 0.86 0.88 0.9 

Sin 2 26 A 

FIG. 2. The figure shows the predictions for sin 2 20A and Atuq/ Am\ for the range of quark 
masses and mixings that fit charged lepton masses. 



model, the Yukawa matrices have a hierarchical pattern, a rough understanding of which 
could come from introducing a local horizontal U(2) H symmetry under which the first two 
families transform as a doublet. This and other aspects of the model such as the inclusion 
of a CP phase are currently under investigation. 

This work is supported by the National Science Foundation Grant No. PHY-0099544. 
We thank K. S. Babu, Z. Chacko and G. Senjanovic for discussions. 
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0.82 0.84 0.86 0.88 

Sin 2 20 A 

FIG. 3. The figure shows the predictions of the model for sin 2 20A and U e % for the allowed range 
of parameters in the model. Note that U e 3 is very close to the upper limit allowed by the existing 
reactor experiments. 
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